The impact and complete genome characterisation of viruses involved in outbreaks of gastroenteritis in a farrow-to-finish holding

Viral enteric pathogens continuously burden intensive pig farming, causing gastrointestinal diseases of epidemic and endemic nature. The present study investigated two diarrhoea outbreaks on a large farrow-to-finish holding and subsequent circulation of outbreak-related enteric viruses. These viruses were characterised by whole genome sequencing, and statistical evaluation of the impact on specific production metrics was performed. The results provided evidence that the Porcine epidemic diarrhoea virus–swine enteric coronavirus (PEDV–SeCoV) S gene recombinant strain was responsible for the first outbreak, whilst Rotavirus A (RVA) in a mixed infection with Rotavirus B (RVB) and porcine kobuvirus (PKV) probably caused the second diarrhoea outbreak. Whole genome characterisation revealed a porcine origin of all viruses involved and significant heterogeneity of RVB strain, proposing four novel genotypes and changes in RVB VP1 genotype classification. The statistical evaluation confirmed only a minor disturbance in the number of weaned pigs per sow, with statistical forecasting showing positive trends. A follow-up study corroborated the endemicity of RVA and PKV, in contrast to PEDV–SeCoV. Punctual, comprehensive and timely investigation of diarrhoea outbreaks is a prerequisite for applying adequate pig health and biosecurity management. Calculating such outbreaks' impact on production metrics can potentially shape future decisions on management improvements.

Rotaviruses (RV) are highly environmentally resistant and ubiquitous, which makes pig health management demanding 15 .These viruses (~ 18.5 kb in size) possess a double-stranded segmented (N = 11) RNA genome and are members of the Reoviride family, genus Rotavirus 3,16 .Out of nine officially recognised species (Rotavirus A-D and Rotavirus F-J) 16 , four are known pathogens in pigs, Rotavirus A (RVA), Rotavirus B (RVB), Rotavirus C (RVC) and Rotavirus H (RVH) 17 .In addition to the standard binomial nomenclature of RV genotypes (VP7 and VP4 segments denominating G and P genotypes, respectively), the whole-genome-based genotype classification system has been established, more precisely Gx-P[x]-Ix-Rx-Cx-Mx-Ax-Nx-Tx-Ex-Hx designate VP7-VP4-VP6-VP1-VP2-VP3-NSP1-NSP2-NSP3-NSP4-NSP5 genomic segments 18 .Rotaviruses in domestic pigs are highly genetically diverse, especially RVA, the most investigated RV species 19 .That was also evidenced in our latest study performed during three consecutive RVA seasons (2018-2021), where we detected eight and seven different G and P genotypes, respectively, in 23 different genotype combinations 20 .In general, rotaviruses in pigs are underinvestigated, especially considering species other than RVA with a lack of whole genome-based studies 17,21 .The clinical outcome of RV infections (indistinguishable from coronavirus infections) ranges from subclinical to severe, resulting in low-to-moderate mortality, mostly in suckling piglets 15 .
Apart from coronaviruses and rotaviruses, many other viruses are involved in the aetiology of neonatal diarrhoea in pigs 22 .For instance, porcine kobuvirus (PKV), classified as Aichivirus C species in the Kobuvirus genus within the Picornaviridae family 23 , is considered a potential causative agent of diarrhoea in pigs 24 , even though its importance as a sole pathogenic agent is debatable 25 .The genome of PKV is a single-stranded positive-sense RNA, only ~ 8.2 kb in size 3 .It circulates in high prevalence in diarrhoeic and healthy pigs 25 .
The present study aimed to investigate the aetiology and impact of two outbreaks of diarrhoea in a large farrow-to-finish holding followed by elevated mortalities in suckling piglets.Subsequent circulation of outbreakrelated enteric viruses was further investigated in a follow-up study.The effect on some production results was measured.Furthermore, the complete genome characterisation of the Porcine epidemic diarrhoea virus, Rotavirus A, Rotavirus B and porcine kobuvirus was revealed.

Sampling
During the first half of 2017, the Croatian Veterinary Institute (HVI) was called to assist in resolving two contagious diarrhoea outbreaks (outbreak 1 in January and outbreak 2 in June) causing elevated mortalities in suckling piglets bred on a large farrow-to-finish holding (1000 sows) in Međimurje County (Northern Croatia).The holding did not implement any vaccination scheme to prevent losses from diarrhoea-causing pathogens.In outbreak 1, diarrhoea, vomiting and mortalities were reported in suckling piglets and diarrhoea, vomiting and agalactia in sows.In outbreak 2, diarrhoea and mortalities were reported in suckling piglets.In January 2017, when outbreak 1 was reported, two faecal samples and three piglet carcasses (age 7-10 days) were sent to the HVI and the Department of Veterinary Pathology at the Faculty of Veterinary Medicine, the University of Zagreb (VEF), respectively.In March 2017, following outbreak 1, three rectal swabs of weaned pigs with diarrhoea were collected and sent to the HVI.In June 2017, when outbreak 2 was reported, four piglet carcasses were sent to the VEF for a necropsy.
The affected holding was further monitored for the circulation of viral pathogens confirmed in both outbreaks in a follow-up study from October 2018 until January 2021.Within that period, 67 rectal swabs were sampled from suckling (N = 55) and weaning (N = 12) pigs and sent to HVI.In total, 79 samples were collected from the affected holding, out of which 84.8% (67/79) originated from pigs with diarrhoea.

Gross and histopathological examination
A routine necropsy examination of seven piglets from two outbreaks (three from outbreak 1 and four from outbreak 2) was performed.Tissue samples (duodenum, jejunum, ileum, colon and mesenteric lymph nodes) collected during necropsy were fixed in 10% neutral buffered formalin, embedded in paraffin, sectioned in five μm thickness, placed on silanised slides and stained with hematoxylin and eosin (HE) for histopathological examination.Intestinal contents were taken from the small intestines of all seven piglets at the necropsy and sent to HVI to perform virological, bacteriological and parasitological tests.

Laboratory diagnostic investigations
Due to the anamnestic data given by the veterinarian working at the farm, the focus was primarily on the molecular diagnostics of PEDV, TGEV and RVA.A 20% w/v intestinal content or rectal swab suspensions were prepared in Medium 199 (Sigma Aldrich, USA), vortexed and centrifuged for 5 min at 14.000 × g.Supernatants were used as starting material for RNA/DNA extraction, which was obtained on KingFisher™ DuoPrime or Flex purification system (Thermo Fisher Scientific, USA) using MagMAX™ CORE Nucleic Acid Purification Kit (Thermo Fisher Scientific, USA) by following the manufacturer's instructions (complex workflow).The exogenous Internal Positive Control (IPC) RNA, Xeno™ RNA Control (Thermo Fisher Scientific, USA), was added to each sample (2 µL) to monitor the appearance of potential PCR inhibitors.The extracted RNA was stored at − 80 °C if not processed immediately.
Detection of PEDV was conducted on two targets, S gene real-time RT-PCR 26 and N gene Sybr Green realtime RT-PCR (unpublished primers, kindly provided by Dr Akbar Dastjerdi, APHA, UK) 6 .TGEV and RVA realtime RT-PCR protocols were focused on detecting the N gene 27 and VP2 genomic segment 28 , respectively.All TaqMan probe assays (PEDV S gene, TGEV and RVA) were performed using VetMAX™-Plus One-Step RT-PCR kit (Thermo Fisher Scientific, USA) and following manufacturer's instructions on reaction set-up and cycling.For all these assays, the final concentration of each primer and TaqMan probe was 600 nmol and 200 nmol, respectively.Due to the dsRNA genome, the real-time RT-PCR set-up for the RVA was divided into denaturation and reaction mixture, as described previously 29 .The QuantiTect Sybr Green RT-PCR kit (Qiagen, Germany) was used for PEDV N gene detection (final primer concentration 500 nmol) in line with the manufacturer's instructions on reaction set-up and cycling.All real-time RT-PCR runs were performed on a Rotor-Gene Q or QIAquant 96 5plex (Qiagen, Germany).
Routine parasitological and bacteriological examinations were used to screen the intestinal contents taken at the necropsy of seven piglets for certain parasites (Isospora suis, Cryptosporidium spp and nematodes) and bacteria (Clostridium spp., E. coli, Salmonella spp.and Lawsonia intracellularis) important from differential diagnostics point of view.Using blood agar and evaluating the presence of haemolysis, the method was able to detect the majority of the enterotoxigenic E. coli (ETEC) F4 and F18 strains 30 .Detection of Lawsonia intracellularis was performed by PCR using GoTaq® G2 Hot Start Colorless Master Mix (Promega, USA) and previously published primers A and B 31 .

Next generation sequencing
The NGS was performed on the MiSeq platform (Illumina, USA) on two samples of small intestinal content.First was the sample from outbreak 1 with the lowest Cq value for PEDV S and N gene real-time RT-PCR, and second was the sample from outbreak 2 with the lowest Cq value for RVA VP2 real-time RT-PCR.Library preparation and the MiSeq run were conducted following previously described procedures 32 .The raw data were analysed using Geneious 8.1.8.software (Biomatters Ltd., New Zealand).After initial data quality control and trimming, virome analysis was performed by Kraken 2 33 and visualised in Krona 34 .Based on virome analysis from the previous step, trimmed reads were mapped to reference genomes of PEDV, RVA, RVB and PKV obtained from the GenBank (PEDV: KU297956, RVA: KX988264-KX988274 and MK936375-MK936425; RVB: KX362400-KX362410 and PKV: KP144318).The mapping was further performed to reference genomes of some other relevant porcine enteric viral pathogens, i.e. astrovirus, Rotavirus C and sapovirus (astrovirus: LC201612, RVC: MT874983-MT874993, MG451776, MG4561777, MG451781 and sapovirus: MK962339).Additionally, de novo assembly was utilised on Geneious 8.1.8.(Biomatters Ltd., New Zealand) software with default settings to assemble all RVA and RVB gene segments due to significant genotype diversity.Consequently, the contigs generated by the de novo assembly were compared for similarity against the whole NCBI GenBank virus nucleotide database using BLASTn.Finally, consensus sequences of PEDV, RVA, RVB, and PKV genomes were analysed for their ORFs using Geneious 8.1.8.software (Biomatters, Ltd., New Zealand), and the deduced amino-acid sequences were obtained.

Follow-up study
All samples taken between 2018 and 2021 (N = 67) and 12 in 2017 during initial outbreaks (intestinal contents from seven carcasses and five faecal/rectal swab samples) were screened for the circulation of the pathogens detected by NGS.We have designed two additional primer sets to detect RVB and PKV based on sequences of corresponding strains determined by NGS in outbreak 2. RVB detection was conducted by primers RVB-VP6-F (5'-TCT GAT CGA GAC AGT GAA TGC-3') and RVB-VP6-R (5'-CTG TGA ACT ACC TGC TCA ATG-3') amplifying 498 bp RT-PCR product of RVB VP6 genomic segment (I9 genotype).PKV detection was performed by primers PKV-3D-F (5'-TGA TTC ACA CTC TGA CAA TG-3') and PKV-3D-R (5'-CGA GAT GTT TCT CAA CAA TG-3') amplifying 507 bp RT-PCR product of PKV 3D gene.Both assays were performed with the utilisation of SuperScript™ III One-Step RT-PCR System with Platinum™ Taq DNA Polymerase (Thermo Fisher Scientific, USA) with reaction set-up as recommended by the manufacturer and the final primer concentration of 600 nmol.Due to the dsRNA genome of the RVB, the reaction mixture was composed of denaturation and reaction mixture.The RNA sample, RVB-VP6-F primer and PCR-grade water were included in the denaturation mixture, followed by incubation at 95 °C for 5 min.After denaturation, a reaction mixture of 2X Reaction Mix, SuperScript III RT/Platinum Taq Mix and RVB-VP6-R primer was added.Both RT-PCR assays (for RVB and PKV) were run on the same cycling conditions on Biometra TRIO (Analytic Jena, Germany) or ABI 9700 GeneAmp thermal cycler (Applied Biosystems, USA) as follows: 30 min at 50 °C; 2 min at 94 °C; 40 cycles of 30 s at 94 °C, 30 s at 50 °C and 30 s at 68 °C; 5 min at 68 °C.
Samples collected between 2018 and 2021 (N = 67) were screened for PEDV (S and N genes), TGEV (N gene), RVB and PKV.RVA detection (VP2 genomic segment) and VP7/VP4 genotyping were performed on those samples within our previous study (19).Since 12 samples directly related to outbreaks were already screened for PEDV, TGEV and RVA, they were additionally checked for RVB and PKV.

Sanger sequencing, genotype assignment, phylogenetic and recombination analysis
The Sanger sequencing was used to confirm the specificity of RVB and PKV RT-PCR amplification.RT-PCR products were purified with ExoSAP-IT™ PCR Product Cleanup Reagent (Thermo Fisher Scientific, USA) as described previously 29 .The samples were subjected to direct Sanger sequencing in both directions using Macrogen Europe (Amsterdam, the Netherlands) services.
Genotype assignment of RVA strain was performed by the BLASTn search in combination with the ViPR tool 35 .Cutoff values for each RVA segment were previously defined 18 .The recently updated genotype classification 36 of RVB strains was combined with the BLASTn to assign a genotype for each RVB segment.
The phylogenetic analysis was done on PEDV, RVA, RVB and PKV strains described in the present study, and a selected number of reference sequences obtained from the GenBank.Multiple sequence alignment was performed by the MUSCLE algorithm, and phylogenetic analysis by the maximum-likelihood (ML) method and models (the lowest BIC score www.nature.com/scientificreports/for VP1-VP3 and NSP1) genomic segment, as well as for the PKV (GTR + G + I).The neighbour-joining (NJ) method combined with the maximum composite likelihood model was used to analyse the PEDV S gene and complete genome.The branching support of the ML and NJ tree was assessed by bootstrap analysis with 1000 repetitions.These analyses were done in MEGA11 software 37 .The phylogenetic trees were visualised and annotated using iTOL 38 .The nucleotide pairwise identity matrix was calculated in BioEdit software, version 7.2.5 39 .
The potential recombination events in the PEDV strain (S gene and complete genome-based) characterised in the present study were analysed using the RDP4 package (version 4.49), which embeds multiple recombination detection methods 40 .Only recombination detected by a minimum of six out of seven methods (RDP, Bootscan, 3Seq, Chimaera, SiScan, MaxChi and Geneconv) was considered valid.The important notice is that "parent" does not necessarily mean the actual progenitor of the recombinant strain but rather a representative of strains whose genome sequences most closely resemble the recombinant 41 .

Immunohistochemistry
Immunohistochemistry (IHC) was used to detect PEDV, RVA and RVB antigens in tissue samples from all seven piglets from outbreaks 1 and 2. Tissue sections were routinely sectioned (5 μm), placed on silanised slides and heated in a heater at 50 °C overnight.After deparaffinisation and hydration, antigen retrieval was performed with citrate buffer (pH 6.0) in an autoclave for 20 min at 121 °C and 1.1 bar.Endogenous peroxidase activity was quenched with 1% hydrogen peroxide in TBST buffer solution for 30 min.Non-specific binding was blocked with normal goat serum (30 min) before incubation with mouse monoclonal primary antibodies to PEDV (SD6-29, Medgene labs, USA; dilution 1:1000), Rotavirus A (3C10, HyTest, Finland; dilution 1:100) and Rotavirus B (10B1, Kerafast, USA; dilution 1:100) for 1 h in a humid chamber at room temperature.PKV monoclonal antibodies were unavailable.A peroxidase-conjugated secondary antibody and detection kit was used for visualisation (REAL EnVision Detection System, Peroxidase/DAB+, Rabbit/Mouse kit, Agilent Dako, USA), after which the slides were counterstained with hematoxylin, dehydrated and coverslipped.

Farm production data
Available time series, i.e., data on piglet mortality, newborns, and weaned pigs per sow, were collected weekly and monthly for three years (2017-2019).Data on the average number of farrowing per sow per year (farrowing index) were also compared.The number of sows for which the calculations were done was 880, 917 and 1017 in 2017, 2018, and 2019, respectively.

Statistics
Statistical analysis was performed to calculate the significance of the impact of two diarrhoea outbreaks on farm production results.Analyses were performed using Systat V.13.2 42 , and subsequent analysis for prediction was performed in Wolfram Mathematica V.13.0.1.0 43.
Infectious disease time series are generally characterised by two separate contributions to the underlying pattern.These components, the trend cycle and seasonal factors can usually be identified by decomposition methods.The trend cycle represents long-term changes, while the seasonal factor denotes the periodic fluctuations of constant length 44 .
The seasonal index (SI) expresses the seasonal component of the time series.If the seasonal index is greater than 1 (or 100%), the incidence is higher than the average and vice versa.The SI (%) was calculated by comparing the average number of newborns, mortalities and weaned pigs per sow weekly in 2017 to the average number in the respective week of all months in 2017.Furthermore, the calculation was expanded to compare these data in each month to the average data of the respective month during three years (2017-2019).After the seasonal indices were calculated, the data were deseasonalised by dividing them by the corresponding indices as follows: deseasonalised data = raw data/seasonal index 45 .Then, the Autoregressive Integrated Moving Average (ARIMA) method 46 was used to estimate the trend based on the 2017-2019 data and to forecast one year ahead (2020).Decisions on statistical relevance were made at the significance level of p ≤ 0.05.

Ethical approval
This study was approved by the Board of Ethics of the Croatian Veterinary Institute (protocol code Z-VI-4-5206/17, approved on 11 December 2017) and complied with ARRIVE guidelines 56 .All applicable international, national, and/or institutional guidelines for the care and use of animals were followed.The pig carcasses were submitted voluntarily by the pig owner and oral consent was obtained prior to sampling.

Laboratory diagnostic investigations
Two faecal samples received in January 2017 tested positive for the PEDV S gene (Cq: 9.9 and 11.7) and negative for the TGEV N gene (SeCoV was excluded) and RVA.Small intestinal contents, taken from three piglets at the necropsy, tested PEDV S gene positive (Cq: 11.8, 17.9, 24.8), PEDV N gene positive (Cq: 9. www.nature.com/scientificreports/Successful detection of IPC RNA excluded the presence of PCR inhibitors.The majority of parasite and bacteria species that are important from differential diagnostic point of view were ruled out as potential causative agents.

Necropsy and histopathology
In all seven piglets, a necropsy showed a small intestine filled with a substantial amount and a large intestine with a smaller amount of watery white-yellowish contents.The intestinal wall, particularly the jejunum and the mesocolon, showed an oedematous appearance.Mesenteric lymph nodes were moderately enlarged.
Jejunum was most severely affected with lesions characterised by shortening, blunting and fusion of villi (Fig. 1).Additionally, in PEDV-infected animals, marked submucosal oedema was noticed (Fig. 1).Duodenum and ileum were affected by similar changes but of lesser intensity.Vacuolisation of superficial enterocytes was mainly observed in PEDV-infected piglets.In RVA-infected animals, the lesions were more intense and additionally characterised by stronger enterocyte necrosis and increased lymphocytes in the lamina propria of the small intestine and colon (Fig. 1).Reactive hyperplasia was visible in the mesenteric lymph nodes.

IHC
Duodenum, jejunum and ileum from all three piglets in outbreak 1 were positive for PEDV antigens.A positive reaction is visible as intracytoplasmic brown granular staining in superficial enterocytes (Fig. 1).Colon and mesenteric lymph nodes were immunohistochemically negative.All samples from four piglets in outbreak 2 were negative for PEDV antigens.All parts of the small intestine (not shown) and the colon from two piglets in outbreak 2 were positive for RVA antigens.Intracytoplasmic brown granular staining was seen in the crypts of Lieberkűhn (Fig. 1).Mesenteric lymph nodes were immunohistochemically negative.All samples from the other two piglets in outbreak 2 and all three in outbreak 1 were negative for RVA antigens.RVB antigens were not immunohistochemically detected.

NGS, phylogenetic and recombination analysis
The NGS data for the sample from outbreak 1 resulted in 566,026 reads, of which 67,561 were mapped to the PEDV reference genome (KU297956).The mean coverage and Q30 were 388.3× and 99.6%, respectively.The NGS data for the sample from outbreak 2 resulted in 568,064 reads, of which 226,467 mapped to referent RVA strains (KX988264-KX988274, MK936375-MK936425), with the mean coverage 1352× and the Q30 of 96.2%.Furthermore, 49,630 reads were mapped to the referent RVB genome (KX362400-KX362410; mean coverage 250×, Q30 95.9%) and 2037 reads to the referent PKV genome (KP144318; mean coverage 30.7×,Q30 95.3%).The www.nature.com/scientificreports/complete CDS were obtained for PEDV, RVA, RVB (except VP2, missing only 2.3%) and PKV.Both de novo and reference mapping assemblies failed to detect the presence of genomes of other viral pathogens in these samples.The PEDV strain responsible for outbreak 1 shared the highest % of nucleotide (nt) and amino acid (aa) identity with the PEDV-SeCoV recombinants reported to circulate in Europe at that time (Fig. 2).The highest identity was with the Polish strain (MZ325485) on the S gene and whole genome level (99.9% on nt and aa levels).The backbone genome originated from the S-INDEL PEDV strain (Fig. 2), with recombination detected only in the S gene by all seven RDP4 methods.The breakpoint positions at 225 (99% CI 208-231) and 636 (99% CI 575-653) of the S gene nucleotide sequence (410 nt long stretch from SeCoV) were identified.The strains KU297956 (SVN/ JH-11/2015) and KX689261 (SVK/42845) were defined as major and minor parents, respectively.The current PEDV strain (HRV 01/2017-Medj) and the Croatian PEDV strain reported in 2016 (CRO/OB-15343/2016) 6 clustered separately (Fig. 2) and shared the S gene nt and aa sequence identity of 98.5% and 96%, respectively.
Compared to the RVA, the lower read numbers were mapped to reference RVB and PKV genomes.For the RVB/Pig-wt/HRV/06-2017-Medj/2017/G27P [6] strain, the complete genotype constellation G27-P[6]-I9-R7-C6-M4-A8-N10-T4-E4-H7 was assembled.Sequence identities to the closest phylogenetic clade representatives are shown in Table 1.The results of the phylogenetic analysis confirm the porcine origin of all genomic segments (Fig. 4, Supplementary Fig. S2).As an important finding, we are proposing the novel genotypes for the VP7 (G27), VP4 (P [6]), VP1 (R7) and VP2 (C6) genomic segments (Fig. 4, Supplementary Fig. S2).While a designation of proposed novel genotypes for VP7, VP4 and VP2 aligns with the current genotype classification 36 , we suggest modification for the VP1 segment.The current RVB strain shared 76.5-79% nt identity with many R4 strains, while the VP1 cutoff is 78% 36 .By increasing the cutoff to 79%, we propose dividing the R4 genotype into three different genotypes (R4 and novel R6 and R7) with excellent bootstrap support (Supplementary Fig. S2).In addition to the evolutionary distinctiveness of the Croatian G27P [6] RVB strain, its heterogeneity was further driven by insertions observed in the NSP4 segment; more precisely, three amino acids (KDK) were inserted at position 194 of the aa sequence.

The impact of diarrhoea outbreaks on farm production results
The results of SI calculation clearly show that the 1st outbreak involving PEDV had only a minor impact on the number of weaned pigs per sow, which was observed in the 4th week of January (− 15.46%) and in the 1st week of February 2017 (− 20.64%) (Fig. 6A).In contrast, outbreak 2, reported in June 2017, had no significant effect on the individual production results (Fig. 6A).On a monthly basis, in 2017, no negative impact was observed (Fig. 6B) for both outbreaks.In subsequent years, no significant changes in the number of weaned pigs per sow were observed (Fig. 6B).Seasonal index calculations showed irregular fluctuations in piglet mortalities and newborns per sow (data not shown).The farrowing index slightly decreased with increasing number of sows: 2.30, 2.27 and 2.20 per sow in 2017, 2018 and 2019, respectively.
Seasonal influence, i.e. irregular periodic fluctuations, was not statistically significant for the number of weaned pigs per sow (p < 0.389).Still, it had an impact on mortalities (p < 0.026) and newborns (p < 0.006).More importantly, the forecasted trend was significantly positive in terms of the number of weaned pigs (p < 0.024) and newborns per sow (p < 0.0001) and favourably negative with respect to piglet mortality per sow (p < 0.0001) (Supplementary Fig. S3).The fit of the ARIMA method was very good, with the residuals showing approximately equal variability along the length of the series and showing no sign of trend or shift.

Discussion
The present study brings new data on the complexity of underlying viral pathogens involved in two diarrhoea outbreaks and their impact on selected production results in a large farrow-to-finish holding.The first outbreak in January 2017 was caused by the PEDV-SeCoV S gene recombinant strain (Fig. 2).Similar PEDV-SeCoV recombinants were involved in diarrhoea outbreaks in Hungary 47 , Slovenia (GenBank), Poland 14 , France (GenBank), The Netherlands (GenBank), Spain 11,13 and Germany 48 (Fig. 2).In addition to conventional molecular methods and NGS, a direct connection to a clinical disease was demonstrated with evident IHC staining of PEDV-SeCoV in enterocytes of affected villi in the small intestine (Fig. 1).Trade of live pigs has been associated with the global transmission of PEDV 49 , which was a probable cause for PEDV-SeCoV introduction into the study holding.The first case of PEDV in Croatia was confirmed a year earlier, in 2016, but that was a typical S-INDEL strain 6 (Fig. 2).Subsequent serological testing in early 2017, revealed the circulation of PEDV in one of the largest pig producers in Croatia which reported significant losses 6,50 .Unfortunately, the sequence of the responsible PEDV strain remained unknown 6 .Although the outbreak reported here coincides in time with that outbreak, it is tempting to speculate that different PEDV variants caused these outbreaks due to significant impact variations.The PEDV-SeCoV recombinant strain described in the present study caused a relatively mild impact outbreak.To the best of our knowledge, our study brings a comprehensive statistical background to this claim for the first time, revealing only a short-lasting disturbance in the number of weaned pigs per sow on a weekly basis (Fig. 6).might have been impacted 9 .Unfortunately, these data were not available for statistical analysis.The farrowing index decreased in three consecutive years whilst the number of sows increased.Possible reasons are diverse, but the restricted capacity of gestation pens largely contributed to reduced welfare, potentially impacting the conception rate.Nevertheless, that was compensated with larger litter size and decreasing mortalities, contributing to positive trends and forecasts.
To conclude, this study provided an insightful view into the background complexity of diarrhoea-causing viral agents in large farrow-to-finish holding and their impact on specific production metrics.Whole genome characterisation of PEDV-SeCoV recombinant, RVA, RVB and PKV strains generated an in-depth look into genome specificities and phylogenetic relationships.The endemicity was corroborated for RVA and PKV, in www.nature.com/scientificreports/contrast to the PEDV-SeCoV recombinant strain.Punctual, comprehensive and timely investigation of diarrhoea outbreaks is a prerequisite for applying adequate pig health and biosecurity management.Calculating the impact such outbreaks had on specific production metrics in the past can potentially shape decisions on management improvements in the future.

Figure 1 .
Figure 1.Histopathological and immunohistochemical findings in PEDV (A-C) and RVA (D-F) infected piglets.Shortening, blunting and fusion of villi in PEDV-infected pigs were most pronounced in the jejunum, accompanied by marked submucosal oedema (A, asterisk) (A).In RVA-infected animals, lesions were also present in the colon, characterised by a moderate increase in the number of lymphocytes in the lamina propria, necrosis of individual surface epithelial cells, and oedema of the mesocolon.(HE 10×) (D).Positive IHC reaction (brown colour) to PEDV and RVA antigens was detected in the intestinal mucosa (jejunum, B, C and colon, E, F) (IHC 10x).Intracytoplasmic brown granular staining was seen in superficial enterocytes in PEDV (C, arrow) (jejunum, IHC 40×) and in crypts of Lieberkűhn of RVA infected animals (F, arrow) (colon, IHC 40x).

Figure 2 .
Figure 2. Phylogenetic relationship between PEDV-SeCoV recombinant strain and selected reference strains.The phylogenetic analysis was based on the S gene (A) and the complete genome (B) nucleotide sequence.The strain described in the present study is bolded and marked in red, while the branch in red designates the cluster of PEDV-SeCoV recombinant strains.The GenBank accession numbers of referent strains are designated within taxa.Both trees were generated by the NJ method and maximum composite likelihood model in MEGA 11 software.The branching stability of each phylogenetic tree was assessed by 1000 bootstrap replicates (values indicated adjacent to the nodes if > 0.7).The scale bar represents the number of substitutions per site.

Figure 3 .
Figure 3. Phylogenetic relationship between RVA strains based on the VP7 (A) and VP4 (B) genomic segments.The strain described in the present study is bolded and marked in red.All sequences are representative of the genotypes G3 (A) and P[23] (B).Human G3 lineages (Lineages I and IX) are annotated.The GenBank accession numbers of referent strains are designated within taxa.Based on the complete CDS, both trees were generated by the ML method and T92 + G + I model in MEGA 11 software.The branching stability of each phylogenetic tree was assessed by 1000 bootstrap replicates (values indicated adjacent to the nodes if > 0.7).The scale bar represents the number of substitutions per site.In displaying RVA strain nomenclature within taxa, the brackets for the P genotype were omitted for the sake of simplicity.

Figure 4 .
Figure 4. Phylogenetic relationship between RVB strains based on the VP7 (A) and VP4 (B) genomic segments.The strain described in the present study is bolded and marked in red, and corresponding novel genotypes G27 and P[6] are annotated in blue.The GenBank accession numbers of referent strains are designated within taxa.Based on the complete CDS, both trees were generated by the ML method and T92 + G + I (VP7) or HKY + G + I (VP4) model in MEGA 11 software.The branching stability of each phylogenetic tree was assessed by 1000 bootstrap replicates (values indicated adjacent to the nodes if > 0.7).The scale bar represents the number of substitutions per site.In displaying RVB strain nomenclature within taxa, the brackets for the P genotype were omitted for the sake of simplicity.

Figure 5 .
Figure 5. Phylogenetic relationship between PKV strains based on the complete genome nucleotide sequence.The strain described in the present study is bolded and marked in red.The GenBank accession numbers of referent strains are designated within taxa.The tree was generated by the ML method and GTR + G + I model in MEGA 11 software.The branching stability of each phylogenetic tree was assessed by 1000 bootstrap replicates (values indicated adjacent to the nodes if > 0.7).The scale bar represents the number of substitutions per site.

Figure 6 .
Figure 6.Seasonal indices (SI) for the average number of weaned pigs per saw in 2017-2019 timeframe.The SI (%) was calculated by comparing the average number of weaned pigs for every saw in each week (A) or month (B) to the average number in the respective week for all months (A) or in the respective month for all three years (B).A significant decrease in the number of weaned pigs per sow is indicated with an asterisk.All months were divided into four weeks, and four additional weeks over the year were added to the nearest week, and the average was calculated.